NOTICE 


When  U.S.  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a  definitely  related  Government  procurement  operation, 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise,  or  in  any  way  licensing  the  holder  or 
any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  be  related  thereto. 


FOREWORD 

Itiis  final  report  was  submitted  by  University  of  Alabama  at  Birmingham, 
Department  of  Chemistry,  Birmingham  AL  on  conpletion  of  contract  F04611-86-K- 
0073  with  the  Astronautics  Laboratory  (AFSC) ,  Edwards  Air  Force  Base,  CA.  AL 
Project  Manager  was  Lt  Pete  Dolan. 

This  report  has  been  reviewed  and  is  approved  for  release  and 
distribution  in  accordance  with  the  distribution  statement  on  the  cover  and  on 
the  DD  Form  1473. 


PETE  DOLAN,  1LT,  USAF 
Project  Manager 


(40^  L. 


jlpdq&Ab— 

SEES# 


STEPE 


L.  RODGERS 


Chief,  Applied  Research  in  Energy 
Storage  Office 


FOR  THE  DIRECTOR 


Division 


'Y  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


'ORT  SECURITY  CLASSIFICATION 

CLASSIFIED 


lb  RESTRICTIVE  MARKINGS 


;urity  CLASSIFICATION  AUTHORITY 


CLASSIFICATION  /  DOWNGRADING  SCHEDULE 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 
Approved  for  Public  Release; 
Distribution  is  unlimited 


:ORMING  ORGANIZATION  REPORT  NUMBER(S) 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AL-TR-90-004 


ME  OF  PERFORMING  ORGANIZATION 
/ersity  of  Alabama 
lirmingham 


6b  OFFICE  SYMBOL 
(If  applicable) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Astronautics  Laboratory  (AFSC) 


DRESS  (City.  State,  and  ZIP  Code) 
irtment  of  Chemistry 
Station 

aingham  AL  35294 _ 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

AL/LSX 

Edwards  AFB  CA  93523-5000 


ME  OF  FUNDING /SPONSORING 
GANIZATION 


8b  OFFICE  SYMBOL 
(If  anr^liratile) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F0461 1-86-K-0073 


DRESS  (City,  State,  and  ZIP  Code) 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

ELEMENT  NO. 

NO. 

NO 

62302F 

5730 

00 

m 

'LE  (Include  Security  Classification) 

High  Energy  Density  Small  Ring  Systems  (U) 


RSONAL  AUTHOR(S) 
nertsma,  K. 


YPE  OF  REPORT 

13b  TIME  COVERED 

14.  DATE  OF  REPORT  (Year.  Month,  Day) 

15.  PAGE  COUNT 

al 

FROM  8609  TO  8911 

9012 

41 

PPLEMENTARY  NOTATION 

>endices  A-F  are  reproduced  by  permission  of  the  American  Chemical  Society  and  Elsevier 
Lence  Publishers,  B.V. 


COSATI  COOES  J 

ELD 

GROUP 

SUB-GROUP 

7 

21 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 
Tetraatomics ,  dicarbides,  diborides,  clusters,  potential 
energy  surface,  energy  storage,  specific  impulse,  bond- 


stretch  isomerism* 

BJreasllaat.’tfc 


B2Li2’ 


B2Be2 ’ 


V 


C2Si2’ 


Li3H, 


HEDM, 


ISTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

k  theoretical  research  investigation  of  four-atom  light  element  species  has  been  under- 
en  in  order  to  find  propellants  with  revolutionary  performance.  The  ab  initio  calcula- 
ts  show  that  the  rhombic  geometry  is  a  global  minima  on  the  potential  energy  surface  of 
y  dicarbides  and  diborides  composed  of  H,  Li,  Be,  B,  C,  and  Si.  By  means  of  bond-stretch 
nerism,  bonding  pattern  change,  these  rhombic  species  can  store  significant  amounts  of 
tgy  -  87  kcal/mole  to  1 22 kcal/mole .  Only  in  the  case  of  B^Be^  was  the  tetrahedral  form 

a  stable.  All  of  the  four-atom  species  found  stable  by  this  investigation  are  estimated 
provide  high  specific  impulses. 

v  r 


5  Accession  Number:  342390 


|Z1  ABSTRACT  SECURITY  CLASSIFICATION  | 

UNCLASS  i  f  IED 

Z2b.  TELEPHONE  (Include  Area  Code) 

22c  OFFICE  SYMBOL 

(805)  275-5760 

AL/LSX 

STRIBUTION /AVAILABILITY  OF  A8STRACT 

UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  Q  OTIC  USERS 


IAME  OF  RFSPON<:ioi  e  i*niv  juai 

i  DOLAN,  2Lt .  USAF 


rm  1473,  JUN  86 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


TABLE  OF  CONTENTS 


INTRODUCTION  1 

Energies  of  Four-Membered  Ring  Structures  1 

Rhombic  Tetracarbon,  C4  2 

Diprotonated  Rhombic  C4  3 

TECHNICAL  DISCUSSION  4 

RHOMBIC  STRUCTURES  -  Structures  and  Energies  4 

Tetraatomic  Dicarbides  4 

Tetraatomic  Diborides  6 

BONDING  PATTERNS  6 

Dicarbides  3 

C4  in  more  detail  11 

Diborides  11 

TETRAHEDRAL  STRUCTURES  13 

ESTIMATED  SPECIFIC  IMPULSES  13 

CLUSTER  MOLECULES  -  Potential  Energy  Surfaces  15 

Disili condi carbide,  Si2C2  15 

Diberylliumdiboride,  B2Be2  17 

Dilithiodiboride,  B2Li2  13 

The  LiaH  Ion-Pair  21 

BOND-STRETCH  ISOMERISM  21 

HYDROGENATED  TETRAATOMICS  23 

Hydrogenated  Diberyllium  Diboride,  B2Be2H2  23 

Dialane(6),  AI2H6,  and  Digallane(4),  Ga2H4  24 

CONCLUSIONS  AND  RECOMMENDATIONS  27 

REFERENCES  29 

APPENDIX  A  -  Publications  A-1 


SUMMARY 


Theoretical  ab  initio  calculations  have  been  carried  out  on  tetraatomic  systems  in  the 
search  for  advanced  chemical  propellants.  In  the  approach  to  identify  systems  with  revolu¬ 
tionary  specific  impulses  (Isp  >  500  sec)  we  concentrated  on  moleculax  systems  that  are  (a) 
composed  of  elements  with  high  o  mbustion  energies,  (b)  that  have  a  high  degree  of  unsat¬ 
uration  (H-deficient),  and  (c)  are  very  strained  (inverted  geometries  of  carbon  and  boron). 
The  research  efforts  were  focussed  on  cyclic  tetraatomics  to  develop  an  understanding  of 
their  structural,  energetic,  and  electronic  properties. 

The  rhombic  forms  of  C4,  C.jH+,  and  C4H22+  were  used  as  parent  molecules.  Although 
the  linear  triplet  C4  form  has  been  observed,  there  are  experimental  indications  for  the 
excistance  of  a  rhombic  isomer  (the  two  structure  are  of  equal  stability).  However,  for  the 
diprotonated  species  the  coexistence  of  the  linear  and  rhombic  forms  in  the  gasphase  is 
unequivocal  and  based  on  charge  separations  that  give  excellent  agreement  with  theoretical 
data.  The  estimated  heat  of  formation  of  726  kcal/mol  for  rhombic  C4H22+  is  one  of  the 
highest  for  a  dication. 

The  question  has  been  addressed  whether  other  elemental  combinations  could  render 
similar  systems  of  high  energy.  Related  to  this  question  we  have  developed  an  under¬ 
standing  of  the  bonding  pattern  of  these  molecules  and  how  their  energy  content  may  be 
manipulated. 

A  survey  of  light  four-atom  systems  at  the  HF/6-31G*  level  of  theory  has  been  con¬ 
ducted  on  dicarbides  and  diborides  with  ligands  ranging  from  Li,  Be,  BeH,  BH,  C,  CH, 
A1H,  to  Si.  Of  the  dicarbides  24  different  rhombic  structures  are  minima  on  the  potential 
energy  surface  and  21  diborides  have  been  characterized  as  are  minima.  While  the  distance 
between  the  bridgehead  carbons  (borons)  is  virtually  constant  throughout  all  dicarbides 
(diborides)  their  bonding  varies  significantly.  To  evaluate  the  viability  of  rhombic  forms 
also  tetrahedral  diboride  combinations  have  been  explored.  Of  the  24  diborides  only  three 
light  elemental  combinations  are  minima.  These  are  B2Li2,  B23e2,  and  B2Be2H-. 

To  inspect  in  more  detail  the  theoretical  levels  required  to  describe  rhombic  tetra¬ 
atomics  and  to  evaluate  their  energic  properties,  full  potential  energy  surface  scans  have 
been  performed  on  C2Si2j  B2Be2,  and  B2Li2.  The  latter  two  light  element  systems  have  es¬ 
timated  Iap  values  of  646  and  531  sec,  respectively.  These  values  arc  based  or.  MF4/6-310* 
binding  energies.  The  calculations  show  that  electron  correlation  has  a  large  influence  on 
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the  relative  stabilities  and  binding  energies,  but  only  marginally  influences  the  geometri¬ 
cal  parameters  of  isomeric  structures.  To  inspect  the  chemical  properties  of  such  systems 
toward  hydrogen,  we  have  investigated  the  potential  energy  surface  of  B2  3e2H2 

Bonding  patterns  have  been  evaluated  by  analysis  of  the  electron  densities.  It  has 
been  found  that  no  transannular  bond  exists  between  the  1.523  Asepaxated  bridgehead 
carbons  of  C4.  However  a  small  reduction  in  length  results  in  CC  bond  formation,  or 
conversely  lengthening  the  CC  distance  results  in  further  depletion  of  electron  density  in 
the  center  of  the  rhombic  structure  as  is  the  case  in  C4H+  and  C4H22+. 

Because  of  the  variation  in  bridgehead  CC  (and  BB)  overlap,  and  the  notable  lack 
of  a  CC  bond  in  C4,  we  have  investigated  whether  breaking  of  the  transannular  CC  (and 
BB)  bond  can  render  vibrational  stable  species  of  higher  energy.  Such  a  mechanism,  bond 
stretch  isomerism,  has  been  found  in  several  systems.  In  the  case  of  C2Si2,  which  is 
isovalent  with  C4,  such  a  rhombic  structure  is  found  78  kcal/mol  higher  in  energy  than 
the  global  rhombic  minimum.  Similarly,  in  the  case  of  B2Be2  a  rhombic  form  is  found 
60.6  kcal/mol  higher  in  energy  than  the  global  minimum.  The  barrier  for  interconversion 
between  the  two  rhombic  structure  has  not  been  determined  but  is  estimated  to  be  high 
because  it  would  represent  a  symmetry  forbidden  process.  This  may  suggest  that  the  high 
energy  rhombic  forms  have  substantial  kinetic  stability. 

Structures  containing  lithium  atoms  like  B2Li2  have  ionic  bonds.  Therefore  the  prop¬ 
erties  of  B2Li2  and  B2Be2  differ  greatly.  In  collaboration  with  Michels  we  investigated 
how  ionic  complexes  may  render  high  energy  materials.  It  was  found  that  in  the  Li3H  the 
H~  anion  is  complexed  to  the  side  of  triangular  Li3+  cation  unit  (i.e.  rendering  rhombic 
Li3H),  thereby  leaving  the  properties  of  this  Li3  +  cation  unit  in  tact. 

An  investigation  exploring  the  potential  of  briding  diboride  units  with  A1  atoms  has 
been  initiated.  Because  of  the  general  lack  of  theoretical  information  available  on  Al- 
A1  interactions,  as  compared  with  that  on  B-B  bonds,  exploratory  studies  have  been 
conducted  on  A12H6  as  well  as  on  the  recently  observed  Ga2H6  and  Ga2H4.  These  studies 
focussed  on  determining  binding  energies  and  the  influences  governing  ionic  versus  covalent 
bonding. 
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INTRODUCTION 


For  decades  chemists  have  been  intrigued  by  small  ring  compounds.  Many  stud¬ 
ies  have  been  devoted  to  highly  strained,  small  hydrocarbons,1  carboranes,2  and  silicon 
derivatives.3  Although  giant  steps  have  been  made  in  materials  science  and  there  is  an 
increasing  understanding  of  Li  clusters,  B  clusters,  Be  clusters,  C  clusters  and  Si  clusters, 
the  information  concerning  bonding  and  energetic  properties  of  mixed  cluster  compounds 
remains  very  limited.  In  the  past  three  years,  we  systematically  investigated  a  senes  of 
over  50  tetraatomic  diborides  and  dicarbides  by  ab  initio  molecular  orbital  theory.  In  this 
report,  we  describe  the  results  of  this  theoretical  survey  study  on  potentially  high  energy 
four-membered  ring  compounds  composed  of  the  elements  Li,  Be,  B,  C,  Al,  and  Si. 

The  foundation  for  developing  high  energy  density  systems  is  often  based  on: 

(a)  increasing  the  unsaturation  by  stripping  molecules  of  all  their  hydrogens. 

(b)  increasing  the  strain  energy  by  inverting  the  carbon  (and  boron)  geometries. 

For  example,  there  are  numerous  strained  hydrocarbons  that  have  either  inverted  tetra- 
coordinated  (sp3)  geometries4  or  pyramidelized  olefinic  bonds.  In  this  investigation,  we 
capitalized  on  the  properties  of  inverted  tricoordinated  (sp2)  boron  and  carbon  atoms.5 
Below  follows  a  description  which  uses  C4  as  an  example  of  how  we  developed  such  un¬ 
saturated  systems  with  inverted  carbon  geometries,  to  find  highly  energetic,  yet  stable, 
strained  ring  systems. 

Energetic  Four-Membered  Ring  Structures;  C4. 

The  expected  increase  in  energy  on  depletion  of  hydrogen  from  the  cyclic  C4 -frame 
is  evident  from  the  HF/6-31G*  data  listed  in  Table  1.  At  this  level  of  theory,  the 
heat  of  formation  A HJ  of  rhombic  C4  (D2/1)  is  228  kcal/mol  (estimated  from  calculated 
heats  of  hydrogenation  in  combination  with  experimental  reference  points — ).  The  deduc¬ 
tion  of  this  value  includes  corrections  for  A  zero  point  energy,  hydrogenation  (each  6.4 
kcal/mol),  A/7/(0-298°),  and  for  differences  (7  kcal/mol)  between  calculated  and  exper¬ 
imental  values.11  The  A HJ  value  of  228  kcal/mol  for  rhombic  C4  (inverted  sp2)  is  large 
compared  to  the  heat  of  formation  of  89  kcal/mol  (HF/6-31G*)  for  stable  [l.l.ljpropel- 
lane  (inverted  sp3),  which  has  a  strain  energy  SE  of  103  kcal/mol.10  The  experimentally 
known  C4H22+  dication,  with  the  highest  A HJ  of  any  cationic  system  (726  kcal/mol),  pro¬ 
vides  further  evidence  of  this  contrast.12  In  this  investigation  we  made  use  of  this  contrast 
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TABLE  1.  6-31G*  energies  'in  a. a.),  heats  of  formation  (AH/,  in  kcal/mol), 


and  strain  energies  (SE,  in  kcal/mol)  for  C4  structures. 


Compound 

Formula 

Energy 

AH/ 

SE 

^4,6 

c4 

-151.14598 

228 

26 

c4h22+ 

-151.53998 

726 

Cyclobutadien«  dication7 

c4h42+ 

-152.91650 

623 

Bicyclo[1.1.0]butene8 

c4h4 

-153.58016 

142 

127 

1 ,3- Cyclobutadiene9 

c4h4 

-153.64120 

103 

65 

Bicyclo[1.1.0]butane8 

c4h6 

-154.87176 

51.9 

65 

Cyclobutene10 

c4h6 

-154.89962 

37.5 

29 

Cyclobutane10 

c4h8 

-156.09703 

6.S 

26 

Hydrogen8 

h2 

-1.12683 

between  C4  and  [l.l.l]propeUane  to  design  over  forty,  stable,  high  energy,  tetraatomic 
species. 

Rhombic  Tetracarbon,  C4. 

We6  and  others13  have  demonstrated  that  rhombic  C4  (1)  is  of  nearly  the  same  energy 
as  triplet  linear  C4.  The  unusual  stability  of  cyclic  C4  can  be  explained  by  inspection  of  its 
molecular  orbital  framework.6  The  highest  occupied  molecular  orbital,  HOMO,  is  a  non¬ 
bonding  cr  orbital  between  the  bridgehead  carbons  with  the  electron  density  located  outside 
the  molecule.  Hence,  the  HOMO  does  not  contribute  to  ring  strain.  The  stabilization  of 
rhombic  C4  results  from  aromatic  7r  bonding;  a  <7-bridged-7r  MO  is  also  present.  Besides  its 
a  skeleton,  the  molecule  has  two  lone  electron  pairs.  In  simple  Lewis  terms,  the  bonding 
of  both  bridgehead  C  atoms  in  C4  can  be  described  as  follows:  (a)  two  valence  electrons 
are  involved  in  ring  C-C  bonding,  (b)  the  third  electron  participates  in  aromatic  ring  n 
bonding,  while  (c)  the  fourth  electron  is  involved  in  a  nonbonding. 
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Scheme  I.  Inverted  tetracoordinate  geometries  at  carbon 


Scheme  II.  Inverted  tricoordinate  geometries  at  carbon 


Diprotonated  Rhombic  C4. 

Our  study  on  the  C4H22  +  dication  underscores  the  validity  of  inverted  planar  tri¬ 
coordinate  carbons.6  Cyclic  C4H22+  (2,  D2<i)  is  only  13.3  kcal/mol  (MP3/6-31G*)  less 
stable  th.  n  its  triplet  linear  diacetylenic  isomer,  despite  the  anticipated  strong  electro¬ 
static  charge-charge  repulsion.12*1  Both  isomers  have  been  observed  in  the  gas  phase. 120 
An  alternative  (but  empirical)  explanation  for  the  stability  of  rhombic  C4H2'+  can  be 
given  along  the  lines  shown  for  strain  build-up  in  saturated  hydrocarbons,  now  adapted 
for  olefins  (Scheme  I  and  II).  Stepwise  replacement  of  the  vicinal  hydrogens  from  the  parent 
ethylene  by  CH+-groups  gives  either  the  cyclopropenyl  cation  or  cyclic  C4H22+,  (note  that 
in  this  series  the  tricoordinate  carbon  is  completely  inverted,  while  the  2-71  stabilization  is 
maintained!). 
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TECHNICAL  DISCUSSION 


The  question  this  investigation  addressed  is  whether  the  inverted  tricoordinate  geom¬ 
etry  concept  developed  for  hi,  h  energy  density  small  ring  hydrocarbons  can  be  expanded 
to  include  other  elements.  There  are  two  important  aspects  to  consider. 

(a)  Do  the  characteristics  of  an  inverted  planar  t.ricoordinate  caxb  jxi  change  using  briuging 
elements  other  than  carbon? 

(b)  Can  the  bridgehead  carbons  be  replaced  by  other  elements  while  maintaining  the  high 
energy  density  properties  of  the  system? 

To  address  these  questions  we  discuss  a  large  number  of  tctraatomic  rhombic  and 
tetrahedral  structures  and  will  analyze  several  small  clusters.  These  latter  potential  energy 
surface  studies  are  presented  in  a  summarized  fashion.  A  complete  z'  fiysis  for  each  cluster 
is  presented  in  the  appendices,  which  have  appeared  or  are  to  appear  in  the  scientific 
literature. 

RHOMBIC  STRUCTURES  —  Structures  and  Energies 

Tetraatomic  Dicarbides.  We  have  examined  at  the  6-31G*  level  of  theory  25 
rhombic  dicaxbides  with  ligands  ranging  Torn  Li,  Be,  BH,  C,  A1H,  and  Si,  in  virtually  all 
combinations.  The  data  are  listed  in  Table  2.  This  group  of  planar  structures  includes 
anions,  cations,  and  dipolar  molecules.  Of  these  25  dicaxbides,  2J,  are  minima  on  the 
6-31G*  potential  energy  surface,  as  indicated  by  the  zero  number  of  imaginary  frequen¬ 
cies.  These  data  strongly  support  the  viability  of  rhombic  structures  that  contain  inverted 
tricoordinate  carbons  irrespective  of  the  type  of  bridging  ligand. 

Similarly  impressive  in  these  dicarbides  is  their  rather  silo'd  C-C  distance  which  has  an 
average  value  of  only  1.457  A  (diamons  excluded).  This  is  unexpected  because  of  the  range 
of  bridging  elements  involved.  Inspection  of  the  Mulliken  overlap  population  indeed  shows 
that  the  nature  of  bonding  between  the  carbons  must  vary  significantly.  Representative 
values  aie  -0.236  for  C4H22+,  -0.063  for  C4,  and  +0.306  for  C2Be2-  However,  the  Mul¬ 
liken  population  analyses  is  subject  to  criticism,  in  particular  for  unsymmetrical  bonding 
arrangements,  and  hence  can  not  rigorously  answer  what  types  of  bonding  are  operative 
in  the  dicarbides.  We  attempted  Weinhold’s  superior  Natural  Bond  Order  approach , 1 4 ■ 1 5 
but  found  too  many  complications  in  using  this  method  due  to  ionicity  in  the  rhombics. 
The  Bader  electron  analysis,  which  we  will  describe  in  the  following  section  on  bonding 
patterns,  allowed  us  to  characterize  the  bonding  in  these  molecules. 
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TABLE  2.  Rhombic  Dicarbides  RHF/6-31G*  Energies,  NIMAG,  and  Geometrical  Data. 
Structure  Bond  Distances 


x-c2-y 

Energy 

Nimag 

c— c 

C— X 

C— Y 

Li-(C2)-Li 

-90.37222 

0 

1.249 

2.108 

2  108 

Li-(C2)-Be 

-97.58752 

0 

1.436 

1.907 

1.571 

Li-(C2)-BH 

-108.30443 

0 

1.426 

1.881 

1.472 

Li-(C2)-CH 

-121.50728 

0 

1.470 

1.968 

1.348 

Li-(C2)-C: 

-120.84757 

0 

1.402 

1.880 

1.418 

Li-(C2)-A1H 

-325.48897 

1 

1.374 

1.922 

1.899 

Be— (C2 )— Be 

-104.79129 

0 

1.459 

1.565 

1.565 

Be-(C2)-BH 

-115.47831 

0 

1.445 

1.558 

1.473 

Be-(C2)-CH 

-128.39388 

0 

1.484 

1.619 

1.345 

Be-(C2)-C: 

-127.98366 

0 

1.418 

1.563 

1.421 

Be-(C2)-A1H 

-332.67417 

0 

1.451 

1.561 

1.857 

BeH-(C2)-BH 

-116.13470 

0 

1.449 

1.670 

1.4.'  ■* 

BeH-(C2)-CH 

-129.26994 

0 

1.501 

1.739 

1.339 

Beh-(C2)-C: 

-128.66892 

0 

1  433 

1.677 

1.408 

BeH-(C2)-AlH 

-333.31024 

0 

1.419 

1.683 

1.853 

BH-(C2)-BH 

-126.15541 

0 

1.470 

1.470 

1.470 

BH-(C2)-CH 

-139.01220 

0 

1.529 

1.526 

1.345 

BH-(C2)-C: 

-138.65329 

0 

1.462 

1.479 

1.418 

BH-(C2)-A1H 

-343.35876 

0 

1.449 

1.470 

1.848 

CH-(C2)-CH 

-151.53998 

0 

1  593 

1.3S6 

1.386 

CH-(C2  )-C: 

-151.45693 

0 

1.520 

1.352 

1.472 

CH-(C2)-A1H 

-356.27607 

0 

1  494 

1.342 

1.919 

C:-(C2)-C: 

-151.14598 

0 

1.457 

1.425 

1.425 

C:-(C2)-A1H 

-355.86570 

0 

1.436 

1.418 

1.857 

A1H-(C2)-A1H 

-560.55536 

0 

1.436 

1.858 

1.858 

:Si-(C2  )-Si: 

-653.41467 

0 

1.415 

1.818 

1.818 

BH-(Si2)-BH 

-628.24903 

0 

2.461 

1.950 

1.950 

:C-(Si2)-C: 

-649.77894 

- 

2.292 

1.863 

1.863 

AlH-(Si2)-AlH 

- 1 057  08620 

- 

2.628 

2.337 

2.337 
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Tetraatomic  Diborides.  We  have  examined  at  the  HF/6-31G*  level  of  theory  28 
rhombic  diborides  with  ligands  ranging  from  Li,  Be,  BH,  G,  A1H,  and  Si,  in  virtually  all 
combinations.  The  6-31G*  data  for  these  systems,  that  contain  two  valence  electrons  less 
than  the  discussed  dicarbides,  are  listed  in  Table  3.  Like  the  planar  dicarbides  this  group 
also  includes  anions,  cations,  and  dipolar  molecules.  Of  the  28  diborides,  21  are  minima 
on  the  HF/6-31G*  potential  energy  surface,  as  indicated  by  the  zero  number  of  imaginary 
frequencies;  the  other  7  are  transition  structures.  Similar  to  the  dicarbides  these  data 
support  the  viability  of  rhombic  structures  now  containing  inverted  tricoordinate  borons 
with  a  large  variety  of  bridging  ligands. 

Considering  the  two-electron  difference  with  dicarbides  these  diborides  have  a  rela¬ 
tively  short  average  B-B  distance  of  1.833  A.  However,  as  will  be  shown,  the  bonding 
nature  of  these  systems  is  entirely  different  from  that  in  the  dicarbides. 

BONDING  PATTERNS 

An  understanding  of  the  bonding  nature  in  the  rhombic  structures  can  be  obtained 
by  analyzing  the  electron  density  ( p )  properties  of  these  molecules.  Prof.  R.F.W.  Bader 
of  McMaster  University  was  so  kind  to  provide  us  with  the  software  that  his  research 
group  has  developed.  A  detailed  description  of  his  method  and  its  applications  can  be 
found  in  the  literature  (also  Appendix  C).16  However,  a  synopsis  of  the  method  assists  in 
comprehending  the  figures  in  this  report. 

In  short,  the  method  is  concerned  with  the  electron  density  p  although  connections 
with  molecular  orbitals  can  be  made.  The  programs  EXTREME  and  PROAIM  were 
interfaced  with  GAUSSIAN  86  and  88  and  use  the  same  basis  sets.  The  properties  of  a 
molecular  charge  distribution,  based  on  Bader’s  topological  analysis  of  density  of  atoms 
in  molecules,  are  summarized  in  terms  of  its  critical  points.  These  are  points  where  the 
charge  density  is  a  maximum,  a  minimum,  or  a  saddle  of  the  gradient  path  of  the  gradient 
vector  field  *.e.,  Vp  =  0.  A  critical  point  is  characterized  by  the  signs  of  its  three  principal 
curvatures  of  p  i.e.,  A].  A2,  and  A3.  A  bond  critical  point  (3,-1)  has  a  minimum  (positive 
curvature)  in  p  on  the  bond  path  connecting  two  nuclei  and  two  negative  curvatures  in 
orthogonal  planes;  the  charge  density  has  the  appearence  of  a  saddle.  Ring  critical  points 
(3,-f- 1 )  have  one  negative  and  two  positive  curvatures.  Three  positive  curvatures  (3, +3) 
define  a  minimum  (cage  critical  point).  The  (3,-3)  critical  point  represents  atoms  or  nuclear 
attractors  ( “pseudo  atoms” )  where  all  curvatures  are  positive.  In  the  figures  presented  the 
connectivity  of  atoms  via  bond  critical  points  are  given  (the  molecular  graphs)  in  both 
plots  of  the  gradient  vector  field,  (which  give  the  separation  of  atoms  in  molecules  by 
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TABLE  3.  Rhombic  Diborides  RHF/6-31G*  Energies,  NIMAG  and  Geometrical  Data 

Structure  Bond  Distances 

X-B2-Y  Energy  Nimag  B— B  B— X  B— Y 


Li-(B2)-Li 

Be-(B2)-Be 

Be-(B2)-BeH 

Be-(B2)-BH 

Be-(B2)-CH 

Be-(B2)-C: 

Be-(B2)-A1H 

Be-(B2)-Si: 

BeH-(B2)-BH 

BeH-(B2)-CH 

BeH-(B2)-C: 

BeH-(B2)-AlH 

BeH-(B2)-Si: 

BH-(B2)-BH 

BH-(B2)-CH 

BH-(B2)-C: 

BH-(B2)-A1H 

BH-(B2)-Si: 

CH-(B2)-CH 

CH-(B2)-C: 

CH-(B2)-Am 

CH-(B2)-Si: 

:C-(B2)-C: 

:C-(B2)-Si: 

:C-(B2)-A1H 

A1H-(B2)-A1H 

A1H-(B2)-Si: 

:Si— (B2 )— Si: 


— 64.04676 
— 7S.3303S 
-78.95048 
-89.03080 
-102.02063 
-101.63581 
-306.22846 
-352.67586 
-89.67334 
-102.86362 
-102.28936 
-306.86504 
-353.32652 
-99.72164 
-112.66014 
-112.28438 
-316.93016 
-363.35948 
-125.29859 
-125.15331 
-329.91335 
-376.30524 
-124.63404 
-375.93214 
-329.53246 
-534.12534 
-580.57306 
-627.00013 


0 

0 

1 

0 

0 

0 

1 

0 

1 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 


1.540 

2.197 

2.197 

2.013 

1.760 

1.760 

1.888 

1.755 

1.862 

1.809 

1.753 

1.598 

1.823 

1.313 

1.443 

1.863 

1.754 

1.383 

1.997 

1  756 

2.040 

1.790 

1.758 

1.911 

1.742 

1.870 

1.594 

1.787 

1.937 

1.445 

1.765 

1.906 

1.39 

1.877 

1.858 

2.036 

1.734 

1.881 

1.920 

1.690 

1.604 

1.604 

1.719 

1.658 

1.444 

1.645 

1.642 

1.399 

1.819 

1.594 

2.032 

1.683 

1.607 

1.941 

1.740 

1.482 

1.482 

1.623 

1.464 

1.460 

1.848 

1.440 

2.102 

1.726 

1.442 

2.047 

1.752 

1.420 

1.420 

1.689 

1.392 

2.004 

1.813 

1.382 

2.037 

1.996 

2.037 

3.553 

1.778 

2.037 

1.910 

1.666 


1.910 


1.947 


“zero  flux”  lines)  and  the  Laplacian  field  of  the  charge  density  V2p  (the  sign  of  which  also 
characterizes  the  polarity  of  a  bond). 

Tetraatomic  Dicarbides.  The  variation  of  bonding  possible  in  planar  rhombic 
dicarbides  can  be  illustrated  in  the  neutral  series  C2Li2,  C2Be2,  C2B2,  and  C4,  which 
shows  a  gradual  changeover  from  the  fully  ionic  C2Li2  system  to  the  covalently  bonded 
C4.  Figure  1  gives  a  representation  of  the  gradient  vector  field  of  their  charge  density 
(HF/6-31G*)  in  terms  of  the  paths  traced  out  by  the  vectors  V/>  and  also  displays  the 
corresponding  contour  maps  of  the  Laplacian  of  their  charge  density  overlaid  with  bond 
paths. 

The  lithiums  in  C2Li2  are  bonded  to  the  center  of  the  C-C  bond.  In  fact,  the  C-C 
bond  path  has  three  bond  critical  points  and  two  “pseudo  atoms”  (nuclear  attractors, 
(3,-3)).  This  unusual  behavior  can  be  interpreted  as  a  C-C  triple  bond  with  a  long  path 
of  virtually  constant  electron  density.  Recently,  pseudo  atoms  were  also  found  in  studies 
by  Bader  and  Gatti  on  Li  and  Na  clusters.17  They  argued  tentatively  that  these  nuclear 
attractors  may  be  associated  with  metallic  behavior.  However,  our  observations  of  these 
(3,-3)  critical  points  in  dicarbides  and  diborides  ( vida  infra )  may  indicate  a  more  universal 
behavior.  In  replacing  Li  for  Be  to  give  CaBej,  bond  paths  are  found  at  the  circumference 
of  the  rhombic  structure,  although  there  remains  a  high  degree  of  ionicity.  The  results  on 
C2Be2  agree  with  those  reported  by  Cremer  and  Schleyer.18  The  difference  between  C2B2 
and  C4  reveals  an  interesting  feature.  While  the  zero-flux  line  of  the  bridging  ligand  in 
C2B2  is  clearly  defining  the  boron,  in  C4  it  is  penetrating  the  rhombic  structure.  In  fact, 
the  two  ring  critical  points  (3,4-1)  and  the  bond  critical  point  (3,-1)  of  the  bridgehead- 
bridgehead  C-C  bond  path  in  C4  nearly  coalesce,  an  observation  also  made  by  Ritchie.19 
Evidently,  the  center  of  C4  represents  a  near  cat  as  trophy  point.  This  means  that  a  small 
perturbation  of  C4  can  eliminate  bonding  between  the  bridgehead  carbons,  despite  their 
short  1.457  A  (HF/6-31G*)  separation. 

The  change  in  bonding  between  the  bridgehead  carbons  becomes  even  more  pro¬ 
nounced  in  comparing  the  series  C4,  C4H+,  and  C4H22+-  Their  6-31G*  transannular  C-C 
distance  are  1.523  (MP2),  1.520  (HF),  and  1.593  A  (HF),  respectively.  The  Laplacian  of 
p  for  these  systems  are  shown  in  Figure  2.  From  this  Figure  it  is  quite  clear  that  there  is 
electron  depletion  in  the  center  of  C4H22+  and  consequently  this  structure  does  not  possess 
a  transannular  C-C  bond. 


8 


Plane 


vector  field  and  the  Laplacian  con 
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lacian  concentration  of  the  HF/6-31G*  charge  density  of  Ca  (MP2).  CaH+,  CdHi2* 


C4  in  more  detail.  Because  of  the  crucial  position  of  C4  in  the  discussed  series  of 
rhombic  dicarbides,  we  wanted  to  study  its  properties  in  more  detail.  As  mentioned,  it 
seems  that  at  the  HF/6-31G*  level  of  theory  the  transannular  C-C  interaction  changes 
from  multiple  C-C  bonding  in  C2Li2  to  no  C-C  bonding  in  C4H22+.  Because  C4  is  the 
structure  where  this  C-  C  bonding  is  lost  (or  formed),  it  seemed  to  us  that  an  analysis  of 
C4  at  higher  levels  of  theory  was  desired.  A  brief  summary  of  this  study  follows.  The  full 
analysis  is  given  in  a  draft  manuscript20  in  Appendix  A. 

The  molecular  orbital  and  topological  electron  density  analysis  of  the  MP2(Full)/6- 
31G*  optimized  rhombic  C4  structure  shows  that  there  is  no  transannular  C-C  bond 
despite  the  1.523  A  separation  between  the  inverted  tricoordinate  carbons.  The  peripheral 
bonds  are  short  (1.451  A,  bond  order  =  1.15)  due  to  strong  surface  (a)  and  7r-delocalization. 
This  contiasts  with  [l.l.l]propellane,  which  at  the  same  level  of  theory,  does  have  bonding 
between  the  inverted  tetracoordinate  carbons  which  are  separated  by  1.592  A. 

Diborides.  To  demonstrate  the  difference  in  bonding  behavior  of  diborides  and 
dicarbides,  we  have  analyzed  the  electron  density  properties  of  the  planar  rhombic  forms 
of  B2Li2,  B2Be2,  and  B4.  Figure  3  displays  the  gradient  vector  field  of  their  charge  density 
in  terms  of  the  paths  traced  out  by  the  vectors  Vp  (left),  while  the  contour  maps  of  the 
Laplacian  of  their  charge  density  overlaid  with  bond  paths  are  given  on  the  right  hand 
side. 

Some  remarkable  features  are  apparent.  In  the  diborides  the  borons  are  bonded  with 
two  bond  critical  points  (3,-1)  located  either  on  the  same  B-B  bond  path  (separated  by  a 
pseudo  atom)  or  on  two  separate  B-B  bond  paths. 

In  the  case  of  the  covalent  B4  structure,  the  four  ring  B-B  bond  paths  have  two 
bond  critical  points  each.  The  center  of  the  molecule  shows  a  depletion  of  charge  with  no 
bonding  between  the  atoms  (instead  it  is  a  (3,+l)  ring  critical  point).  The  bonding  pattern 
in  planar  B2Be2  is  different.  The  two  borons  are  bonded  by  two  curved  bond  paths  with 
an  “electron  density  hole”  (i.e.,  a  (3,-1- 1)  ring  critical  point)  in  the  center.  The  molecule 
contains  four  conventional,  although  slightly  polarized  B-Be  bonds.  A  combination  of 
some  of  the  B4  and  B2Be2  features  is  found  in  B2Li2.  This  cluster  has  two  highly  curved 
B-B  bond  paths,  each  containing  two  bond  critical  points.  The  center  between  the  borons 
also  contains  a  bond  critical  point  which  connects  the  “pseudo  atoms”  on  the  two  B-B 
bond  paths,  to  which  also  the  two  lithiums  are  bonded.  This  structure  has  no  bonding  at 
the  circumference  and  is  highly  ionic. 
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Molecular  Graphs  of  Diborides  at  the  Molecular  Plane 
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This  analysis  and  the  comparison  with  diborides  makes  clear  that  these  tetraatomic 
diboride  clusters  have  special  bonding  features,  which  were  as  of  yet  undiscovered.  In  our 
opinion  we  have  just  scratched  the  surface  of  the  bonding  features  possible  for  boron  based 
clusters.  A  more  detailed  study  of  B2Be2  and  B2IU2  is  given  in  section  2.5.2  and  2.5.3, 
respectively. 

TETRAHEDRAL  STRUCTURES 

Light  elemental  clusters,  like  the  tetraatomic  Be4,  may  form  tetrahedral  structures 
( vida  infra).  Because  several  rhombic  species  were  calculated  to  be  transition  structures 
we  decided  to  evaluate  tetrahedral  diboride  structures,  which  can  also  be  considered  as  high 
density  species.  An  incentive  for  this  study  was  also  provided  by  the  fact  that  the  above 
discussed  rhombic  structures  axe  not  necessarily  the  global  minima.  For  example,  linear 
triplets  can  be  energetically  competitive  as  in  C4,  C4H+,  and  C4H22+,  (all  16  valence 
electrons)6  and  even  be  strongly  preferred  as  in  the  12  valence  electron  C2Be2-18  The 
lighter  10  valence  electron  ionic  dilithiodicarbide  C2Li2,  however,  has  again  a  rhombic 
structure.21  Although,  Li4  is  also  a  planar  rhombic  (in  sharp  contrast  to  H4),  a  higher 
energy  tetrahedral  form  is  also  feasible. 

We  calculated  24  tetrahedral-like  diborides  at  the  HF /6-31G*  level  of  theory.  Their 
geometries  and  energies  are  listed  in  Table  4.  Of  the  24  structures  S  are  minima  (B2Li2, 
B2Be2,  and  B2Be2H~)  and  16  axe  transition  structures.  This  data  shows  that  densily 
packed  tetrahedral  structures  are  preferred  for  certain  light  elemental  combinations. 

ESTIMATED  SPECIFIC  IMPULSES 

Having  established  that  rhombic  carbon  structures  axe  of  high  energy,  and  having 
determined  that  indeed  a  large  variety  of  tetraatomic  rhombic  and  some  tetrahedral  clusters 
are  minima  (but  not  necessarily  global  minima)  on  the  potential  energy  surface,  a  crude 
estimate  of  their  specific  impulses  can  be  obtained  with  the  following  formula: 

Isp  =  265  x  >/[(  heat  of  formation)/ (mol.  weight)] 

The  heat  of  formation  is  estimated  from  our  ab  initio  calculated  binding  energy  and 
the  heats  of  formation  of  gaseous  atoms  from  elements  in  their  standard  states.  In  Table  5 
ab  initio  energies  (MP4  and  HF/6-31G*)  of  rhombic  structures  are  used  except  for  B2Be2 
for  which  the  more  stable  tetrahedral-like  isomer  was  used  (vida  infra). 
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TABLE  4.  Tetrahedral  Diborides  RHF/6-31G*  Energies,  NIMAG  and  Geometrical  Data. 


Structure 

x-b2-y 

Energy 

Nimag 

B— B 

Bond  Distances 

B— X  B— Y 

X— Y 

Li-(B2)-Li 

-64.06043 

0 

1.548 

2.252 

2.252 

3.504 

Be-(B2)-Be 

-78.37354 

0 

1.548 

1.772 

1.772 

2.156 

Be-(B2)-BeH 

-79.01726 

0 

1.529 

1.802 

1.926 

2.447 

Be-(B2)-BH 

-88.99520 

1 

1.520 

1.764 

1.653 

2.131 

Be-(B2  )-CH 

-101.87523 

1 

1.447 

1.962 

1.717 

1.571 

Be-(B2)-C: 

-101.50896 

1 

1.484 

2.102 

1.511 

1.557 

Be-(B2)-A1H 

-306.26094 

1 

1.544 

1.766 

2.049 

2.554 

Be-(B2)-Si: 

-352.64423 

1 

1.524 

1.763 

1.993 

2.569 

BeH-(B2)-BH 

-89.64323 

2 

1.531 

1.903 

1.625 

2.226 

BeH-(B2)-CH 

-102.74791 

2 

1.525 

2.173 

1.4S0 

1.852 

BeH-(B2)-AlH 

-306.90465 

1 

1.535 

1.927 

1.994 

2.752 

BH-(B2)-BH 

-99.59929 

2 

1.567 

1.656 

1.656 

1.908 

BH-(B2)-CH 

-112.51562 

1 

1.427 

1.884 

1.675 

1.450 

BH-(B2)-A1H 

-316.87552 

1 

1.528 

1.660 

2.038 

2.404 

BH-(B2)-Si: 

-363.24008 

2 

1.445 

1.844 

2.259 

1.945 

CH-(B2)-CH 

-125.16082 

1 

1.447 

1.716 

1.716 

1.340 

CH-(B2)-C: 

-125.00759 

1 

1.430 

1.682 

1.687 

1.342 

CH-(B2)-A1H 

-329.75801 

1 

1.500 

1.502 

2.246 

2.129 

CH-(B2)-Si: 

-376.15976 

1 

1.433 

1.654 

2.332 

1.826 

:C-(B2  )-C: 

-124.61891 

2 

1.407 

1.750 

1.750 

1.383 

:C-(B2)-A1H 

-329.40172 

1 

1.488 

1.516 

2.405 

1.814 

A1H-(B2)-A1H 

-534.14341 

1 

1.538 

2.043 

2.043 

2.912 

A1H-(B2  )-Si: 

-580.52553 

1 

1.526 

2.040 

1.997 

2.823 

:Si— (B2 )— Si: 

-626.90122 

1 

1.458 

2.135 

2.135 

2.364 

TABLE  5.  Estimated  Isp  values.0 
System  Binding  Energy  Specific  Impulse 


(kcal/mol) 

MP4  HF 

( 

MP4 

sec) 

HF 

B2Li2 

178.5 

87.8 

591 

727 

B2Be2 

195.4 

122.7 

646 

73S 

c4 

576.4 

496.8 

399 

524 

C2Si2 

379.0 

415.5 

396 

354 

These  data  confirm  our  supposition  that  rhombic  structures  indeed  have  large  esti¬ 
mated  Iap  values  and  that,  as  expected,  they  are  largest  for  those  systems  that  are  com¬ 
posed  of  elements  with  the  highest  combustion  energies.  The  data  for  B2Li2  and  B2Be2 
also  demonstrate  the  significant  influence  of  electron  correlation  in  determining  calculated 
binding  energies.  Therefore,  the  given  l3p  values  are  only  rough  estimates,  because  of 
the  direct  relationship  between  the  estimated  magnitude  of  the  specific  impulse  and  the 
accuracy  of  the  calculated  binding  energy. 

CLUSTER  MOLECULES  —  Potential  Energy  Surfaces 

In  this  section  we  elaborate  in  more  detail  the  potential  energy  surfaces  for  a  selected 
group  of  tetraatomics.  Although  we  have  demonstrated  that  most  rhombic  structure  are 
viable  species  (minima),  it  is  also  apparent  that  they  need  not  be  the  global  minima.  To 
further  evaluate  the  potential  of  rhombic  (or  tetrahedral)  species  for  consideration  as  high 
energy  density  systems,  it  is  relavant  to  investigate  more  fully  the  structures  of  certain 
elemental  combinations.  We  report  here  our  analysis  of  C2Si2,  B2Be2,  B2Li2,  and  L^H, 
which  cover  a  group  of  widely  different  species  ranging  from  light  to  heavy  elemental 
combinations. 

Disilicondicarbide,  Si2C2.  Our  first  objective  was  to  investigate  a  binary  cluster 
for  its  structural,  energetic,  and  electronic  properties.  As  already  noted,  the  16  valence 
electron  C4  shows  a  slight  preference  for  the  rhombic  over  linear  structure,  which  reverses 
on  mono-  and  diprotonation  to  C4H+  and  to  C4H22+.6  In  contrast,  the  related  16  elec¬ 
tron  Si4  strongly  favors  a  rhombic  form,22  with  its  20  electron  tetraanion  SU4-  fragment 
being  tetrahedral.23  It  is  well  recognized  that  the  bonding  properties  of  silylenes  differ 
significantly  from  those  of  the  corresponding  carbon  compounds.3  It  this  then  relavant  to 
investigate  the  properties  of  a  mixed  cluster  that  contains  both  C  and  Si.  We  give  here 
a  very  brief  summary  of  our  investigation  of  the  potential  energy  surface  for  the  C2Si2 


system.  A  full  description  of  this  work  has  appeared  in  the  literature2'1  and  is  attached  in 
appendix  B. 

The  structures  and  energies  of  the  binary  disilicon  dicaxbide  C2Si2  in  the  lowest  singlet 
(s)  and  triplet  (t)  states  (14  isomers)  were  investigated  by  ab  initio  MO  theory.  Full 
fourth-order  Moller-Plesset  (MP4)  perturbation  theory  was  employed  on  seven  HF /6-31G* 
optimized  geometries.  Rhombic  dicarbide  3s  is  the  global  minimum  and  8.5  kcal/mol 
more  stable  than  the  rhomboidal  structure  4s,  which  contains  an  inverted  tricoordinate 
carbon  as  well  as  an  inverted  tricoordinate  silicon.  Both  cyclic  structures  are  energetically 
favored  over  a  linear  triplet  isomer  5t.  Whereas  the  relative  energies  are  very  sensitive 
to  electron  correlation  effects,  addition  of  diffuse  functions  (MP2  and  HF/6-31-+-G*)  has 
liotle  influence. 


1  723 


1  268 


Si 


C - C 


Si 


5t 

This  study  shows  that  disilicon  dicaxbide  has  similarities  with  the  related  C4  and  Si4 
but  also  that  it  is  inherently  more  complex;  an  important  aspect  of  this  will  be  emphasized 
in  section  2.6.  We  note  that  there  is  an  isoelectronic  relationship  between  3s  and  the 
rhombic  A^S^2-  fragment  in  the  CaA^Sb  material,25  which  is  representive  for  a  large 
class  of  solid  Zintl  complexes.26  The  rhomboidal  structure  4s,  with  its  mixed  C/Si  bonding, 
is  also  unique  to  binary  clusters  and  similarly  may  well  be  an  important  structural  type. 
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Diberylliumdiboride,  B2Be2-  This  10  valence  electron  species  is  composed  ot  the 
two  elements  with  the  highest  combustion  energies.  Because  of  its  anticipated  importance 
as  a  high  energy  material,  B2Be2  was  subjected  to  a  full  (singlet  s  and  triplet  t)  potential 
energy  surface  study  to  characterize  its  global  minimum,  to  determine  its  stability,  and 
to  investigate  its  binding  energies  and  properties.  Because  of  its  smaller  difference  in 
elemental  electron  negativities  than  that  in  the  isoelectronic  dilithium  dicarbide  (C2Li2  has 
a  preferred  (6°)  puckered  rhombic  structure)21  and  the  12  val.  el.  diberyllium  dicarbide 
C2Be2,18  the  properties  of  B2Be2  were  expected  to  deviate  from  C2IU2  C2Be2,  and  other 
tetraatomics.  Beryllium  borides  are  known  since  the  turn  of  the  century  with  crystal 
structures  (powder  patterns)  reported  for  BBeg,  BBc3.  BBe2,  B2Be,  B4Be,  BgBe,  BgBe, 
BioBe,  and  B^Be  with  the  noticable  exception  of  B2Be2-2‘  Here  we  give  here  a  very  brief 
summary  of  our  study  of  the  B2Be2  system.  A  full  analysis  is  given  in  Appendix  C  and  is 
to  appear  in  the  literature.28 

The  ab  initio  molecular  orbital  study  of  the  10  val.  el.  dibervllium  diboride  (B2Be2) 
potential  energy  hypersurface  yields  a  singlet  tetrahedrallike  form  6s  as  the  global  mini¬ 
mum  at  the  HF  and  MP2/6-31G*  levels  of  theory.  At  MP4/6-31G*  the  energy  difference 
of  6s  with  any  other  isomer  is  at  least  30  kcal/mol;  (high  spin-contamination  in  triplet 
structures  was  amhilated).  The  topological  electron  density  analysis,  based  on  the  theory 
of  atoms  in  molecules,  shows  6s  to  have  its  Be  atoms  polar  tt-  complexed  to  the  3-B  bond 
with  an  angle  of  80.1°  between  the  ~  planes;  despite  the  short  2.156  A  Be-Be  distance, 
there  is  no  Be-Be  bond.  The  inversion  barrier  (6s  ^  7  s  ^  6s)  amounts  to  11.7  kcal/mol 
at  MP4/6-31G*  +■  ZPE  and  14  6  kcal/mol  at  MP2/6-31  lG(df)  +ZPE.  The  D2/1  st  uctures 
7s  and  8s  are  related  by  means  of  “bond-stretch’"  isomerism.  The  electron  density  analysis 
shows  that  the  electron  density  is  concentrated  in  the  -'enter  of  structures  7s  whereas  8s 
shows  a  depletion  of  electron  density  in  the  center  as  compared  to  the  periphery  of  this 
rhombic  structure.  The  trapezoidal  9s  and  the  Be-substituted  B2Be  ring  10s  are  disc  rssed 
for  their  structural  and  electronic  properties.  All  B2Be2  structures  show  short  B-B  bonds 
(t.e.,  1.549  A  for  6s,  MP2/5-31G*)  with  the  exception  of  8s  (1.869  A,  MP2/6-31G*).  The 
effects  of  electron  correlation  on  relative  energies  and  binding  energies  is  significant.  The 
binding  energy  of  the  four  elements  in  6s  s  204  kcal/mol  and  that  of  one  beryllium  in  6s 
is  76  kcal/mol  (all  MP4/6-31G*).  A  calculated  specific  impuls  for  6s  is  given  in  section 
2.4. 
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MI‘2/6  31G* 
HF/6-31G* 
(11F/3-21G) 


7s 


I.S4S 
(1  556) 


B 


1.785 
(1  800) 


Be _ Be 

1  996 
(2.015) 


9s 


Be 


10s 


From  this  study  we  infer  that  the  nature  of  bonding  (covalent  and  ionic  complexation) 
is  an  important  aspect  for  the  formation  of  laxge  clusters  that  may  form  molecular  solids. 
This  supports  the  finding  of  Dr.  Brener’s  at  the  Louisiane  State  University  that  larger 
clusters  (molecular  solids)  of  the  tetrahedral  B2Be2  unit  cell  are  favorable. 2tl 

Dilithiodiboride,  BaLU-  Both  finely  dispersed  Li  and  boron  compounds  are  known 
to  boost  the  performance  of  rocket  fuels.  Therefore,  our  next  step  was  to  investigate  the 
properties  of  ‘he  even  lighter  tetraatomic  system,  the  8  val.  el  dilithio  diboride  BtLU- 

An  ab  initio  molecular  orbital  study  was  performed  on  the  B2LU  singlet  and  triplet 
potential  energy  hypersurface.  The  singlet  rhombic  form  11s  is  the  global  minimum  at  the 
HF  and  MP4/6-31+G*  levels  of  theory  of  the  10  isomers  investigated.  Figure  4  displays 
the  HF/6-31G*  optimized  structures.  At  the  correlated  level  this  isomer  is  favored  by  8.3 
Kcal/mol  over  the  tetrahedral  isomer  and  by  at  least  20  kcal/mol  more  than  any  other 
isomer.  However,  the  data  listed  in  Table  6  show  a  rather  strong  oscillation  between 
the  MP2,  MP3,  and  MP4  relative  energies  for  11s  and  12s,  indicating  that  probably  a 
rnultireference  wave  function  is  required  for  determination  of  accurate  energies.  The  effects 
of  electron  correlation  on  all  relative  energies  is  very  significant  as  graphically  shown  in 
Figure  5.  A  calculated  specific  impuls  for  11s  is  given  in  section  2.4. 
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TABLE  6.  Total  (au)  and  Relative  (kcal/molc)  Eneigics  of  Li2B2  Isomers 


Structure 

HF/6  31G* 

Rel.  En.  NIMAG  <s2> 

MP4/6-31G* 

Rel.  En. 

11s 

-64.04676 

0.0 

-64  30754 

0.0 

lit 

-64.05044 

-2.3 

(0)  2.03 

-64.25693 

31  8 

12s 

-64.06043 

— S.6 

;o 

-64.29335 

8.9 

12t 

-64.08114 

-21.6 

1)  2.25 

-64.26819 

24.7 

13s 

-64.05724 

-6.6 

°) 

-64.24349 

40.3 

13t 

-64.03463 

7.6 

2)  2.03 

-64.23644 

44.6 

14s 

-64.01299 

21.2 

>2) 

-64.24051 

42.1 

15t 

-64.07447 

-17.4 

0)  2.21 

-64.25147 

35.2 

16t 

— 64.0S134 

-21.7 

[0)  2.05 

-64,26682 

25.6 

17t 

-64.06390 

-10. S 

1)  2.05 

-64.25268 

34.4 

<•4.116  (4.032)»- 


11 


2.252^^* 

(2-275) 


^1548 

6/(1580) 


8  a504 

(2JB63) 


12 


1591  2.135 

Li — Bfissfe-fifissUi 

13 


Bond  lengths  for  triplet  species 
ire  given  in  parenthesis. 


Figure  4.  HF/6-31G*  geometries  of  B9L19  isomers. 
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The  topological  electron  density  analysis  for  the  rhombic  B2Li2  structure  shows  that 
the  Li  atoms  complex  ionicly  to  the  B2-unit.  This  is  evident  from  the  gradient  vector  field 
and  the  Laplacian  of  p,  which  are  displayed  in  Figure  6.  The  molecular  graph  for  11s, 
also  displayed  in  Figure  6,  shows  a  complex  topology.  This  propably  results  from  the  use 
of  a  limited  basis  set  (MP2/6-31+G*)  to  describe  the  electron  density  properties  of  the 
electron-rich  B22-  group.  The  bond  properties  for  11s  and  12s  are  listed  in  Table  7. 

The  LiaH  Ion-Pair  This  study  was  performed  in  collaboration  with  J.A.  Mont¬ 
gomery  and  H.H.  Michels  of  United  Technology  Reseach  Center.  Our  interest  in  this  4  val. 
el.  system  relates  to  the  chemistry  of  small  clusters  and  ion-pair  states.  A  full  analysis  of 
our  study,  which  has  appeared  in  the  literature,30  is  given  in  Appendix  D.  A  brief  abstract 
of  the  work  follows. 

Ab  initio  calculations  on  LiaH  ion-pair  states  and  a  topological  analysis  of  the  charge 
density  of  the  resulting  optimized  structures  are  reported.  The  global  LiaH  minimum  is  a 
planar  C2V  structure  with  H-Li3+  ion-pair  character.  The  pyramidal  Ca„  structure,  which 
is  21.5  kcal/mol  higher  in  energy  (MP2/6-311++G**),  is  found  to  have  ionic  H-Li  bonds, 
but  no  Li— Li  bonds.  Although  vibrational  analysis  indicates  this  structure  to  be  a  true 
minimum  on  the  potential  energy  surface,  inclusion  of  the  vibrational  zero-point  energy 
makes  the  thermodynamic  stability  of  this  structure  questionable. 

BOND-STRETCH  ISOMERISM 

In  the  course  of  this  work  we  determined  that  many  tetraatomic  systems  may  have 
the  properties  of  bond-stretch  isomerism.  We  resorted  to  the  C2Si2  and  B2Be2  systems 
to  analyze  this  phenomenon  in  more  detail  and  to  establish  its  value  in  small  cluster 
molecules.  The  properties  of  the  B2Be2  system  have  been  fully  analyzed  in  the  already 
discussed  B2Be2  potential  energy  surface  study  that  is  given  in  Appendix  C.28  The  study 
of  the  bond  stretch  isomerism  in  C2Si2  is  discussed  separately  in  full  detail  in  Appendix 
E.  A  brief  abstract  of  this  published  research31  follows. 

The  differences  in  bonding  and  electronic  properties  between  the  MP2/6-31G*  opti¬ 
mized  disilicon  dicarbide  structures  3  and  18  can  be  related  to  bond  stretch  isomerism, 
which  occurs  through  a  level-crossing  mechanism.  Their  energy  differences  axe  calculated 
to  be  78.4  at  the  MP4  and  91.4  kcal/mol  at  GVB/6-31G*.  Because  of  the  symmetry  for¬ 
bidden  transformation  no  barrier  for  this  process  was  calculated.  Electron  density  analysis 
indicates  concentration  of  charge  in  the  center  of  rhombic  structure  3  with  some  multiple 
bonding  character  between  its  transannular  carbons.  In  contrast,  the  higher  energy  isomer 
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Figure  6.  Molecular  graph,  gradient  vector  field,  and  the  Laplacian  concentration  of  the 
MP2/6-31+G*  charge  density  of  B9L1-  ills'). 


TABLE  7.  MP2/6-31+G*  Bond  Properties  for  11s  and  12s 


Struct. /Crit.Pt. 

Type 

Ai 

A2 

A3 

e 

dir 

p(r) 

VV(r) 

H(  r) 

11s  a 

[3.-1] 

-1.21 

-0.40 

0.32 

2.09 

1 

0.73 

-1.28 

-0.60 

b 

[3-3] 

-0.67 

-0.64 

-0.58 

0.04 

II 

0.73 

-1.89 

-0.57 

c 

[3,-1] 

-1.90 

-0.65 

0.67 

1.94 

1 

0.73 

-1.88 

-0.61 

d 

[3,-1] 

-0.80 

-0.43 

4.66 

0.87 

II 

0.21 

3.43 

0.02 

e 

[3, +3] 

0.28 

1.10 

5.33 

— 

- 

0.71 

6.70 

-0.56 

f 

[3,+l] 

-1.55 

0.76 

3.20 

— 

— 

0.71 

3.41 

-0.59 

12s  a 

[3,-3] 

-5.87 

-3.10 

-0.33 

0.89 

1 

1.09 

-9.29 

-0.95 

b 

[3,-1] 

-6.35 

-3.05 

2.92 

1.08 

JL 

1.08 

-6.4S 

-1.13 

c 

[3,-1] 

-0.65 

-0.23 

3.63 

1.80 

-L 

0.18 

2.75 

0.01 
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18  has  a  depletion  of  charge  in  the  center  relative  to  the  circumference  of  the  rhombic 
structure. 

HYDROGENATED  TETRAATOMICS 

Another  question  we  addressed  concerns  the  effect  of  hydrogen  substitution  on  small 
inorganic  clusters.  In  contrast  to  the  extensive  literature  on  carbon  systems,  very  little  is 
known  about  the  related  inorganic  systems.  In  light  of  the  special  bonding  characteristics 
associated  with  B2Be2  (section  2.5.2)  we  decided  to  investigate  the  effect  of  its  hydrogena¬ 
tion.  Since  boranes  are  energetic  materials  and  because  little  is  known  about  alanes  we 
also  explored  Al2Hg.  The  latter  research  is  a  first  step  to  a  closer  investigation  of  rhombic 
and  tetrahedral  diborides  with  bridging  aluminum  ligands. 

Hydrogenated  Diberyllium  Diboride,  B2Be2H2.  An  ab  initio  molecular  orbital 
study  was  performed  on  the  B2Be2H2  singlet  potential  energy  hypersurface.  Structure  19 
is  the  global  minimum  of  the  eight  isomers  analyzed.  These  eight  HF/6-31G*  optimized 
B2Be2H2  structures  are  shown  under  Table  8  that  lists  their  total  and  relative  energies 
at  the  HF  and  MP4/6-31G*  levels.  At  the  correlated  level  19  is  energetically  favored  by 
at  least  14  kcal/mol  over  any  other  isomer.  It  is  also  evident  that  the  effect  of  electron 
correlation  is  much  less  pronounced  than  was  the  case  in  the  hydrogen- void  B2Be2  system. 
However,  the  geometries  and  the  data  listed  in  Table  8  also  indicate  special  properties 
related  to  the  hydrogenation  of  tetrahedral/rhombic  B2Be2.  For  example,  the  energy 
difference  between  hydrogenation  on  the  borons  versus  hydrogenation  on  beryllium  (which 
formally  become  tricoordinate)  amounts  to  only  15  kcal/mol. 

It  appears  that  hydrogenation  of  B2Be2  is  preferred  on  both  borons  followed  by  hydro¬ 
genation  on  the  two  berylliums.  The  global  minimum  19  can  be  regarded  as  a  diborane(2) 
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complexed  with  two  beryllium  atoms.  Consequently,  structure  19  would  also  be  a  very 
energetic  system.  We  anticipate  that  similar  to  the  tetrahedral-like  B2Be2  isomer  also  19 
will  not  contain  Be-Be  bonding  (despite  the  1.886  A  seperation)  but  instead  will  have  its 
berylliums  complexed  to  the  center  of  the  1.764  A  B-B  bond.  The  ‘second  best’  structure 
it  ca.  15  kcal/mol  less  stable  than  19  and  contains  a  1.5  A  short  B-B  bond  with  bridging; 
Be-H  groups  in  either  a  planar  (22)  or  ‘butterfly’  (21)  orientation;  the  nature  of  the  sta¬ 
tionary  point  must  be  determined  at  correlated  levels.  Clearly,  the  bonding  characteristics 
of  this  system  warrant  further  investigation. 

Dialane(6),  AI2H6,  and  Digallane(4),  Ga2H4.  In  contrast  to  the  vast  amount 
of  literature  on  small  boranes,  carbon,  and  silicon  molecules,  only  very  little  is  known  on 
small  aluminum  based  molecules.  Since  aluminum,  for  a  second  row  element,  has  a  rather 
high  combustion  energy  and  because  of  its  likely  suitability  (as  boron)  to  form  high  density 
materials,  we  decided  to  investigate  the  experimentally  still  elusive  dialane(6).  An  addi¬ 
tional  incentive  for  this  study  was  the  recent  experimental  observation  of  digallane(6).32 
Therefore,  a  theoretical  study  on  AI2H6  and  Ga2H6  would  enable  us  to  calibrate  the  ab 
initio  data  against  experimental  Ga2H6  values  as  well  as  to  determine  the  binding  energy 
and  bonding  properties  of  A12H6.  Since  we  recently  reported  on  dialane(4)  A12H4,33  it  was 
also  decided  to  study  Ga2H4  with  the  intention  to  further  calibrate  hydrogenated  group 
III  elements.  This  research  is  considered  to  be  a  necessary  first  step  for  closer  inspection 
of  rhombic  and  tetrahedral  systems  that  contain  aluminum  ligands.  A  brief  summary  of 
our  studies  on  Al2He/Ga2H634  and  Ga2H435  follows.  The  full  analysis  of  these  systems  is 
given  in  Appendix  F  and  G,  respectively,  which  are  to  appear  in  the  literature. 

The  structure,  harmonic  frequencies,  and  binding  energies  of  dialane(6)  and  digal- 
lane(6)  are  calculated  with  ab  initio  MO  theory.  Excellent  agreement  is  found  with  IR 
spectroscopic  data  on  Ga2Hg.  The  binding  energy  for  AI2H6  is  31.5  kcal/mol  at  MP4/6- 
31G**  and  for  Ga2H6  a  value  of  23.4  kcal/mol  is  obtained  with  a  Huzinaga  valence  triple 
(  basis  set  after  correction  for  ZPE  and  electron  correlation  effects  at  MP4. 

The  potential  energy  surface  of  digallane(4)  Ga2H4  has  been  studied  by  ab  initio 
molecular  orbital  theory  at  the  MP2  level  using  a  Huzinaga  valence  triply  (  basis  set  sup¬ 
plemented  with  d-polarization  functions.  Six  of  the  eight  HF/3-21G*  optimized  structures 
were  characterized  as  minima.  Their  bonding  properties  were  studied  by  a  topological 
electron  density  analysis.  Of  the  ionic,  covalent,  and  /j-hydrido  bridged  species  the  tri- 
and  bidentate  Ga+[GaH4]-  forms  27  and  28,  respectively,  are  the  most  stable  isomers. 
However,  the  energy  differences  with  the  other  isomers  are  less  than  15  kcal/mol.  The 
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TABLE  8.  Total  (au)  and  Relative  (kcal/mole)  Energies  of  B2Be2H2  Isomers. 


Structure 

Tot.  En. 

19 

c2u 

-79.62424 

20 

D2/! 

-79.56060 

21 

c2t, 

-79.59930 

22 

d2A 

-79.58958 

23 

c3 

-79.58598 

24 

c3 

-79.55500 

25 

c3 

-79.53958 

26 

c2v 

-79.53812 

Rel.  En.  NIMAG 


0.0 

(0) 

39.9 

(0) 

15.7 

(0) 

21.8 

(1) 

24.0 

(0) 

43.5 

(1) 

53.0 

(0) 

54.0 

(1) 

Tot.  En. 

Rel.  En 

-79.93443 

0.0 

-79.89878 

22.4 

-79.91033 

15.1 

-79.91144 

14.4 

-79.85155 

52.0 

-79.84456 

56.4 

-79.83565 

62.0 

-79.83335 

63.4 

H  132-5 

1-83' 


—  B 


H 

'1.186 


^-^---1886 

Be 


B^x'-B  =  68.4 


19  C2„ 


/N 

B«— X— B«  =  101.2 


20  D2h 


23  C, 


22  D2fc 


X-8-8 


123. * 


24  C3 


25  C. 


26  C2u 


25 


H-..„  / 
.Gai - Gaj 


H 


H 


\ 


H, 


29  D2a 


v  Gai - Hj - Gaj 

H '  i 

H 


30  O3W 


\Ga2 - Gaj - H< 

h-7 

H 

31  C3v 


The  potential  energy  surface  of  digallane(4)  Ga2H4  has  been  studied  by  ab  initio 
molecular  orbital  theory  at  the  MP2  level  using  a  Huzinaga  valence  triply  (  basis  set  sup¬ 
plemented  with  d- polarization  functions.  Six  of  the  eight  HF/3-21G*  optimized  structures 
were  characterized  as  minima.  Their  bonding  properties  were  studied  by  a  topological 
electron  density  analysis.  Of  the  ionic,  covalent,  and  ^/-hydrido  bridged  species  the  tri- 
and  bidentate  Ga+[GaH4]~  forms  27  and  28,  respectively,  are  the  most  stable  isomers. 
However,  the  energy  differences  with  the  other  isomers  are  less  than  15  kcal/mol.  The 
covalent  D2j  isomer  29  is  3.5  kcal/mol  less  stable  than  27.  The  head  and  tail  complexes 
of  :GaH  with  GaH3  (30  and  31),  which  are  of  similar  energy,  are  ca.  13  kcal/mol  higher 
in  energy  than  29.  The  potential  energy  surface  of  Ga2H4  is  similar  to  that  of  A12H4  but 
differences  are  noted. 


26 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  main  conclusions  from  our  general  survey  on  small  cluster  molecules  (ca  60)  and 
from  our  studies  on  seven  selected  systems  are  as  follows: 

The  ab  initio  molecular  orbital  study  on  tetraatomic  systems  has  shown  that  24  of 
the  25  investigated  rhombic  dicarbides  and  21  of  28  studied  rhombic  diborides  axe  local 
or  global  minima  on  the  potential  energy  structures.  This  confirms  the  general  viability 
of  rhombic  clusters  composed  of  the  elements  H,  Li,  Be,  B,  C,  Al,  and  Si.  We  found 
that  only  the  rhombic,  tetrahedral,  and  (triplet)  linear  structures  compete  for  the  global 
minima.  Only  3  tetrahedral  (light)  diboride  clusters  out  of  24  combinations  axe  equilibrium 
structures.  We  have  also  demonstrated  the  high  energy  content  in  these  rhombic  and 
tetrahedral  clusters,  which  axe  composed  of  high  combustion  light  elements.  The  most 
energetic  systems  investigated  (MP4/6-  3lG*//MP2/6-3lG*)  axe  B2Li2  and  B2Be2-  They 
have  the  highest  estimated  Iap  values  (591  seconds  and  646  seconds,  resp.). 

The  general  survey  study  also  establishes  that  the  bonding  character  between  the 
bridgehead  atoms  (B  or  C)  in  rhombic  and  tetrahedral  structures  varies  widely  (t.e.,  from 
double  bonding  to  no  bonding)  and  depends  on  the  nature  (electronegativity)  of  the  bridg¬ 
ing  ligands.  For  example,  rhombic  C4  has  no  transannular  bonding,  but  C2Be2  has  strong 
C-C  bonding  as  does  C2Si2-  Similarly  B4H2  has  no  transannular  bonding,  but  in  B2L12 
there  is  a  multiple  B-B  bond. 

The  potential  energy  surface  studies  on  B2Li2,  B2Be2,  L^H,  and  C2Si2  demonstrate 
an  extreme  sensitivity  of  these  tetraatomics  to  the  effects  of  electron  correlation. 

Bond-stretch  isomerism  is  a  likely  mechanism  for  energy  storage  and  is  common  to  the 
frame  of  rhombic  tetraatomics  as  demonstrated  in  B2Be2  and  C2Si2-  Such  isomers  have 
large  energy  differences  and  axe  separated  by  level  crossing  mechanisms,  implying  that  the 
high  energy  isomers  lie  in  deep  potential  wells. 

To  study  the  structural,  energetic,  and  electronic  properties  of  small  clusters,  com¬ 
posed  of  light  elements,  the  total  electron  density  of  such  systems  must  be  analyzed  in 
combination  with  standard  ab  initio  molecular  orbital  analysis. 
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The  main  recommendations  from  this  research  toward  high  energy  density  materials 

are: 

(a) .  Higher  levels  of  theory  are  required  as  well  as  a  mapping  of  all  dissocation  modes  for 

the  light  elemental  clusters  (H,  Li,  B,  and  Be)  to  obtain  kinetic  and  thermodynamic 
data  of  enough  accuracy  to  predict  reliable  Iip  values. 

(b) .  More  work  is  required  to  predict  the  stabilities  and  properties  of  small  clusters  in 

convential  fuel  systems. 

(c) .  An  in-dmth  study  of  the  scope  of  bond-stretch  isomerism  is  needed  both  to  evaluate 

its  potential  as  an  energy  promoting  mechanism  and  to  find  practical  systems  that 
serve  as  energy  storage  devices. 
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